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We report on a study of the quasielectron-quasihole and skyrmion-antiskyrmion 
bound states in the v = 1 quantum Hall regime. The short range attraction po- 
tential is assumed to be determined by a point magnetic impurity. The calculations 
are performed within the strong field approximation when the binding energy and the 
characteristic electron-electron interaction energy are smaller than the Landau level 
spacing. The Excitonic Representation technique is used in that case. 

I. Unique properties of the two-dimensional electron gas (2DEG) in strong magnetic fields 
attract much attention to its spectrum. In particular, the interaction of 2DEG with neutral 
short range impurities exhibits unusual features in comparison with its 3D prototype.^ In 
this paper we study the bound fermion states appearing in the quantum Hall "ferromagnet" 
(QHF) regime; i.e. the filling factor is v — M /N^ ~ 2n+l, where Af and = L 2 /2tt1 2 3 are 
the numbers of electrons and magnetic flux quanta (L 2 is the 2DEG area, Ib is the magnetic 
length). In the high magnetic field limit, which really represents the solution to the first 
order in the ratio r c = (e 2 / 'els) /huj c considered to be small (u c is the cyclotron frequency, e 
is the dielectric constant), we get the ground state with zeroth, first, second,... and (n— l)-th 
Landau levels (LLs) fully occupied and with nth level filled only by spin-up electrons aligned 
along B. 

In the clean limit, fermion excitations are classified by their spin-numbers \AS Z \ = K+l/2 
(K is an integer),— ranging from the simplest |A5^| = 1/2 case of quasielectrons or quasiholes 
to the K^-oo limit which corresponds to the so-called skyrmions. Certainly, the total energy 
of excitations incorporates the Zeeman energy le^AS^, and the spin number of lowest-lying 
fermions is thus determined by the actual value of the gap ez = \gnsB\- As regards to the 
Coulomb exchange energy of the fermions, in the v = 1 case this part of the total energy 
decreases monotonically with the K number. 3 E.g., to the first order in r c the exchange 
energies are — \\Jir/2e 2 /slB and |y 'n /2e 2 j 'sIb (in the strict 2D limit) for electron-like and 
hole-like skyrmions, respectively.— i 5 i 6 i 7 i 8 For comparison, in the K = case these ones are 



(electrons) and a/vt /2e 2 /els (holes). 

We calculate the bound-state energies in the presence of point impurities. In the single 
electron approximation (i.e. in the v<^ 1 case) this problem was studied in Ref. y for arbi- 
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trary strengh of the perpendcular magnetic field. In Refs. the authors investigated 

the genesis of the impurity potential in the two-dimensional channel and gave a proof for the 
short range approach. The latter in terms of the "envelope" wave-function method means 
that the inpurity Hamiltonian can be modeled by the 5-function: 

Amp = 2vr \Wi + jdz^j 5(t) (1) 

(here and in the following cr x ,y,z are the Pauli matrices). The D ^ case corresponds to a 
paramagnetic impurity with its own magnetic moment aligned parallel to the magnetic field. 
(We consider that B || z.) At the negative g-factor the positive value of D provides capture 
of the spin wave (spin exciton) by an impurity in the QHF case. 9 

We will solve the problem in the "shallow" impurity approximation, 

\W a \«h 2 /m: (a=U), (2) 

where W^i = W±D/2. Actually one will see that this condition enables to employ the high 
magnetic field approach: 

E b < hu c , (3) 

where is the desired binding energy, and u c is the cyclotron frequency. The condition 
(3) allows us to employ the projection onto a single LL approach when calculating E b in 
the leading approximation. At the same time we still consider that the point- impurity 
approximation is not disturbed, i.e. the impurity localization radius pb is much smaller than 
the magnetic length: p& <C 

II. Any single-electron state may be presented in the form of the expansion 

X = C ^a V -, (4) 

ap 

where we choose the Landau-gauge functions 4> ap as the basis set. The subscript p distin- 
guishes between different states belonging to a continuously degenerate Landau level, and 
the label a is a binary index a = (n a ,a a ), which represents both LL index and spin index. 
We have thus <p ap (r,a) = 5^ aa (l B L)' 1/2 e ipy ip na (plB + x/l B ), where 



{2 n n\^)- 1 ' 2 e- x2l2 H n {x) (5) 



[H n (x) is the Hermite polynomial]. In the following we employ the notation a p , b p ,... for the 
electron annihilation operator corresponding to sub level a, b,... and use also the intra-LL 
"displacement" operators A^, <B+ where 

K = N f l E e ~^H + ? a p-f ' - (6) 
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(Aq = A-q). Considering the quantity x as an annihilation operator in the Schrodinger 
representation, we can substitute Eq. (4) into (x|-ffi mp |x) to obtain the secondary quantized 
representation of the Hamiltonian Tii mp , namely: 

n imp « i B 2 e ~ q2ll/A i w ^X + w iK) ■ (7) 

q 

Here the sign of approximate equality means that we have omitted the terms corresponding 
to the LL mixing and have kept only those relevant to the projection onto the nth LL. 
Therefore, specifically we have in Eq. (7) that the labels a and b correspond to a = (n, t) 
and b = (n, j) sublevels. At u = 2n+l the "clean" ground state |0) is completely determined 
by the equations -4. q |0) =5 q ,o|0) and £> q |0) = 0, while electron and hole are defined as the 

l/e)=EWl°) ^ lA> = E/*(P) a pl°> ^ 
V V 

states, respectively. Here the envelope functions are normalized as \f e ,h\ 2 = 1- 

The quasiparticle states (8) satisfy the "clean" equations Ho\fe,h) = (Eq + E et h)\f e ,h), 
where E e = ez and Eh = Eq- Here Eq is the characteristic Coulomb energy: Eq = 
J e -9 2/ l / 2 v (q)dq, where 27rl / (q) is the 2D Fourier component of the averaged Coulomb 



potential (in the ideal 2D case V = e 2 / nq and Eq — -Jit /2e 2 ' j 'kIb)- To obtain this result, it 
is convenient to employ the expression for the Coulomb interaction Hamiltonian in terms of 
the Excitonic Representation (see, e.g., Ref. 0) and the commutation rules 

op] = —e-^^m^ [B+, 6+] ^ j-e-^ B ^l%^ y (9) 

and [A+,b+] = [B+,a p ] =0. 

First we obtain the correction to the ground state energy in the case of single impurity. 
Substituting expression (7) for Hi mp into the equation (Ho + fii mp j |0) = (E + AE )\0) 
(where 7io is the "clean" QHF Hamiltonian including the Zeeman and Coulomb interaction 
part) we obtain the correction to the "clean" ground state: AE = W^/l 2 B . 

Finally, to calculate the bound electron state we solve the equation 7^i mp |/ e ) = E' e \f e ), 
which with help of Eqs. (7)-(9) is reduced to the integral equation 

% / dsf e (s)e-^ l2 B ' 2 = E'J(p). (10) 



The latter has the solution f e (p) = (ttJV^-VV"^/* at E' e = WJ1 2 B [(vr^)" 1 / 4 is the 
normalization factor]. Hence the electron binding energy is 

e { : ] = -wjii. (ii) 
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Naturally, this state is realized if Wj_ < 0. Let us note that in the leading approximation 
the obtained E^ e ' value is equal to the binding energy in the single electron problem.— In 
a similar way we obtain that the equation Hi mp \fh) = E' h \fh) has the solution fh(p) = 
(nN^^e-P 21 ^/ 2 at E' h = -W^/l 2 B . So, the total energy of the |/ h ) bound state is E +E c ; 
i.e. it is just the same as in the "clean" hole state. The bound energy of the hole is 

4 h) = wM- (12) 

This state exists under the condition W-\ > 0. 

The physical meaning of the envelope functions f e ^ obtained above becomes evident if 
we change the Landau gauge to the symmetric gauge. In the latter case we have to change 
the vector potential A= (0, Bx, 0) for A— (—By/2, Bx/2, 0). Then the single electron states 
of the n's LL are described by the basis spatial function 

1/2 

{ir/l B ) m L™{r 2 /l 2 B 2)e- im *- r2/Al B (n + m > 0), (13) 

where r = (r cosy?, r sin tp), L™ is the Laguerre polynomial and m runs over integer 
numbers: m = —n, 1 — n, 2 — n, N^ — n — 1. All these states have the same cyclotron 
energy, and the Fermi creation and annihilation operators acquire index m (instead of p in 
the Landau gauge) now. For example, one can find the expression for old Landau gauge 
operator a p in terms of new operators a m : 

M-l 

a p = Nf 1 ' 2 ^ l Vm(pl B )a m -n (14) 

m=0 

where <p m is the oscillatory function (5). Now, if we substitute this expression (and analogous 
one for 6+) into Eqs. (8) with the above found functions, we obtain that only the m = 
harmonic of Eq. (13) contributes to the \f e ,h) bound states. Indeed, the summation over p 
in the expansion (8) turns out to be proportional to the integral e~ p2 H m (p). The latter 
vanishes at any m except for m = 0. Certainly, this feature reflects the well known fact that 
the point impurity Hamiltonian is diagonal exactly in the symmetric basis (13). Besides, 
only the zero (axially symmetric) harmonic contributes to the bound state energy calculated 
within the single LL approximation. 

III. To study a bound skyrmionic excitation we present it in accordance with Refs. QJsl 
as a smooth rotation in the three-dimensional spin space 

^(r) = U (r)x(r), r = (x, y) . (15) 



nl 



2 m + 1 (m + n)\iT 
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Here ip is a spinor given in the stationary coordinate system and \ is a new spinor in the 
local coordinate system following this rotation. The rotation matrix U(r) (U^U = 1) is 
parameterized by three Eulerian angles:— 

* ( cos f e -^ + ")/ 2 sin § e ^-v)/ 2 \ 

U= 2 2 . (16) 

^ - sin f e^)/ 2 cosfe^")/ 2 y 

These angles 0(r), ip(r) and 7/(r) present continuum field functions. The skyrmion state is 
thus determined by the continuum matrix U(r) and by the local quantum state x determined 
in terms of small gradient ( /^V?7) corrections to the local QHF, where all electron spins 
are parallel [x oc („)]. 

The solution may be found from the reformulated variational principle. Namely, we divide 
the 2DEG area into a large number of domains which are much smaller than the total 
2DEG area but still remain much larger than the magnetic flux quantum area 2itl 2 B . The 
energy of excitations of this type may be found through the minimization procedure in the 
following way: 



E = min 

u 



mm 



(17) 



Here averaging is performed over the domain Gi. All the Gi areas add up to the total 2DEG 
area. Hi is the Hamiltonian corresponding to the Gi domain. The state presents here 
a many-electron quantum state built by the single electron spinors (15). So, the state \ip) 
is parametrized by U and by the derivatives of U (generally, up to any order) considered as 
external parameters for every Gi. The procedure of the inner minimization in Eq. (17) is 
equivalent to the solution of the Schrodinger equation within the area Ax Ay = Gi. 

In our case the local Hamiltonians are just the same as in the clean case except for the 
Hamiltonian H i=Q corresponding to the one domain Gq which involves the point impurity. 
The procedure of the minimization (17) only differs from the clean case by adding the energy 

AEimp = — m^, — (18) 

to the sum within the square brackets in Eq. (17). Here we have used Eq. (15) and then 



should substitute 



U+H imp U = 2 . 2 S(t) (19) 

-£sin0e^ W-^cos9 



2 ^"""^ " 2 



into Eq. (18). (We have set rj = (p without any loss of generality.) Under the conditions 
(2)-(3), it is sufficient to use as \x) the unperturbed QHF state determined in the local 
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coordinate system of domain G . That is, in our case \x) = x|0), where x is the annihilation 
operator. If employing again the expansion of Eq. (4), then with the help of equation 



H imp = J drx U+H imp Ux (20) 

over Go 

we can obtain the impurity Hamiltonian in terms of the secondary quantized representation. 
We will study only the v — 1 case. Substituting Eqs. (19) and (4) into Eq. (20) and 
considering that a = (0, f) and b = (0, 1), we find within the single LL approximation 

^im P = l B 2 E q e-" 2 ' |/4 [{W + f cos 6 )A+ + (W - § cos 6 )B+ 
-f sin# e-^iV; 1/2 g q - § sin^'iV; 1/2 Q+] . 

Definitions of the operators and £>+ are identical to those given by Eq. (6) with the 
number of magnetic flux quanta = Ax Ay /2nl 2 3 being non-zero only for the domain 
G . The notation of the spin-exciton creation Q+ = N^ 1 ^ 2 J2 P e ~ l9xpl2B b^ + qy a p -^. an d 
annihilation Q q = (Qq) + operators have also been used in Eq. (21). 9 and </? are the 
Hermitian angles (determined by the given U matrix) corresponding to the domain Go- Now 
in accordance with Eq. (18) we get 

AE^l = (0\H imp \0) ={w + ^ cos^o) /ll- (22) 

(Only the ~ Aq term contributes to this result.) 

The outer minimization in Eq. (17) is thereby presented as min (^Eu + AE^^j, where 
Ejj is the "clean" energy obtained for a given function U (r) after the summation over all Gi 
in Eq. (17). Meanwhile the minimization of E v and AE^ may be fulfilled independently. 
Indeed, the clean skyrmion state is degenerate having the energy which does not depend on 
the "skyrmion center" position. The latter is the point of the total 2DEG area where 9 = it, 
i.e. local electron spins are aligned in the direction opposite to the direction at the infinity 
(where = 0). In contrast to this, the energy AE^ just depends on the relative positions 
of the impurity site (the r=0 point in our coordinate system) and the skyrmion center site. 
Thus the impurity lifts the degeneracy. Hence by varying the position of the skyrmion center 



site, we find that min 

u 



AE iu) 

imp 



{W-\D\/2)/l%. If D < 0, this value is reached at 6 = 0, 
i.e. the skyrmion center is located at the infinity. Evidently in this case the skyrmion does 
not form a bound state. At D > 0, the minimum energy is realized for 9 = n, i.e. where 
the impurity site and the skyrmion center coincide (of course, to within the length smaller 
than the characteristic skyrmion radius R* but perhaps larger than /#). This means that the 
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bound skyrmion state takes place. The binding energy should be obtained by subtraction of 
the minimum energy from the 60 = energy (W+D/2)l^, i.e. 

Ef = D/l%. (23) 

So, in the adopted approximation only the magnetic impurity with D > captures the 
skyrmion (cf. Ref. |9| where the similar result is obtained in the case of the bound spin 
exciton). It is worth to note that the result (23) does not depend on the skyrmion charge: 12 
electron-like and hole-like skyrmions have the same bound energy (23). 

Analysis reveals that the charge dependence arises only in the second order approximation 
in terms of Z_bV, being determined by the relative correction (Ib/R*) 2 - This can be easily 
found by means of the renormalization procedure for the magnetic length Ib — > Ib- Indeed, 
the effective local magnetic length is determined by the effective magnetic field, including 
the additional part proportional to the second-order spatial derivatives of the field U:& 



L B L B 



V x Q z . (24) 



Here 

-(l + cos0) V<p (25) 

(see Refs. ffl aad I. is the « length at the ^ (far ho m the skyrmlo n eeate r) . 
Let us assume that the impurity is non-magnetic, i.e. D = 0. The binding energy (if the 
bound state would be realized) should be determined exactly by the desired correction. We 
substitute into Eqs. (24)-(25) for the angles their expressions in terms of functions of r — 

R* 2 - r 2 

cos8 = — 5 , and (p = — q arctan(y/a;) 

R* + r 2 

(where q = ± is the skyrmion charge). Eventually, by comparing the skyrmion energy at 
r = with the energy at the infinity, we find the impurity correction: AE imp = WV x 
Q z = 2qW/R* 2 . The binding energy is thus 

r=0 

Ef = -2qW/R* 2 if D = . 

Therefore the bound electron/hole-like skyrmion arises when W is positive/negative. 

In reality, the concentration of point impurities in the 2D channel can be considerable. 
At least it seems to have the values at which the mean distance between impurities is well 
shorter than the effective skyrmion radius (the latter is determined by small but still non- 
zero Zeeman gap e* z ). In this case the impurity contribution to the total energy is equal to 
J dr\lg 2 [W + ^ cos#(r)] , where A(r) is the concentration of impurities. It involves also the 
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correction to the ground state energy l B 2 (W+^ ) J dr\ which should be subtracted. The 
impurity correction to the proper skyrmion energy is thereby 

A£it P = ^- / [cos0(r)-l]Adr. 

B J 

If we compare this value with the skyrmion Zeeman energy 

^ J (l-cos0)dr/2nl 2 B , 

it becomes evident that the magnetic impurities dispersed in the 2D channel lead to the 
correction to the Zeeman gap: 

e* z -> e z - 2ty\D 

(assuming a homogeneous concentration A). Due to the supposed smallness of e* z , this can 
be substantial. Under certain conditions it could change the sign of effective (/-factor in the 
2DEG. 
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